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O ' Abstract 

o : 

The Fermilab crystal-extraction experiment E853 at Tevatron was simulated 
by Monte Carlo code CATCH [1] tested earlier in the CERN-SPS experiment 
[2-4] . Predictions for the extraction efficiency, angular scans and extracted beam 
profiles are presented. Several ideas are proposed and tested by the simulation, 
, how to get in E853 the key information of the extraction experiment: the "septum 

Q I width" of a crystal and dependence of extraction efficiency on it, the impact 

' parameters of protons at crystal, and the contribution of the first and multi 

passes to the extraction. The ways to optimize E853 are analyzed. 
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■ 1 Introduction 



The crystal-extraction experiments at CERN SPS [2-4] and Fermilab Tevatron 
[5] have in view possible application of channeling for proton extraction from 
a multi- TeV machine [6]. The technique employs a bent crystal placed in the 
beam halo, which traps and bends the particles parallel to the crystallographic 
plane within Lindhard angle Oc- The halo particles hit a crystal very close to its 
edge, with impact parameter h in the range ~A to ~/im. This calls for a good 
perfection of the crystal edge. Alternatively, one should investigate how crystal 
extracts particles in the multipass mode, which involves several scatterings in the 
crystal of the circulating particles. 

As the extraction includes many passes, there is no easy way to extrapolate 
the results with energy. This makes the detailed comparison of the measurements 
with computer simulation essential. Such an analysis made [7] for the CERN-SPS 
experiment has shown good agreement of the theory with measurements [2,3]. 
The major result of [7] was a prediction of the edge imperfection of the crystals 
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used at SPS. The new SPS experiment, employing a crystal with an amorphous 
edge-layer to testify this idea, has proved much the same efficiency indeed [4]. 
Another prediction, made for the "U-shaped" crystal - much the same efficiency 
but narrow (70 ^urad fwhm) angular scan[7] - has also been confirmed [3]. With 
the simulation code [1] tested at SPS, here we model the extraction of 900-GeV 
protons from Tevatron, matching E853 [5,8]. 

The real crystal has an irregularity of the surface, which defines a range of 
inefficient b at the edge (" septum width" t) . The following information is essential 
for understanding the crystal extraction process: 

(a) efficiency F, and contributions to it from the first/multi passes; (b) distri- 
bution over h at the crystal; (c) septum width i; (d) dependence F{t). 

We propose the ways to get this information in E853. 



2 Qualitative discussion 

The essential feature of E853 is the fact that the crystal atomic planes are per- 
pendicular to the crystal face in touch with the beam. In E853 one should align a 
crystal in 2 planes: the channeling plane (vertical, y) with the accuracy of Oc, and 
the horizontal plane (x) to keep the crystal face parallel to the incident protons 
(Fig. 1). At first glance, the need to tune two angles is an inconvenience. Here 
we show that this extra degree of freedom is an excellent possibility to study 
extraction in many details! 
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Figure 1: Horizontally tilted crystal: (a) aligned, (b) tilt x' >0, (c) x' <0. 

Consider a crystal disaligned from the beam horizontally at x' . Depending on 
the sign of x', cither the upstream end approaches the beam (we define x' <0), 
or the downstream one {x' >0); fig. 1. Because of x' , a septum width as thick 
as t = \x'\L occurs at the crystal edge; L is the crystal length (4 cm in E853). 
Protons incident in the range 0<6< t do not traverse the full length of crystal. 
The result depends dramatically on x' sign. 

In case of x' > 0, protons traverse the downstream edge. It is disaligned by 
~0.64 mrad (the bending angle) w.r.t. the beam. Therefore, protons traverse it 
like an amorphous substance. This case imitates a crystal with an amorphous 
near-surface layer as wide as t ~ x'L. Measuring F{x') for x'>0, one measures 
F(t). Theory [9] predicts very weak F{t) dependence at high energies. The 
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confirmation would be encouraging for tfie multi-TeV crystal extraction. Notice 
that the step of t scan could be very fine: with Sx'=2.5 /xrad and L=4 cm one 
has St= 0.1 ^m. 

In case of x' < 0, protons traverse the upstream edge which is aligned w.r.t. 
the beam. Therefore, many particles are trapped in channeling. However, those 
incident in the range 0<b<x'L traverse a reduced (<4 cm) length, thus getting 
a reduced (<0.64 mrad) deflection and therefore are lost. The inequality of two 
cases, x'>0 and x'<0, causes a strong asymmetry of F{x') dependence. The 
difference /S.F=F{x')~F{—x') is proportional to the number of protons incident 
with 0<6<x'L. Varying x' and observing AF, one investigates the distribution 
over h at crystal, with accuracy of (56=0.1 jum. 

This is complicated by another interesting phenomenon. The protons incident 
on an aligned imperfect crystal with hmax < t, have to traverse the full length 
of the crystal, and to experience a substantial nuclear scattering. Suppose, this 
crystal is disaligned so that bmax/x'L rjO.1. Then at first incidence the protons 
traverse only the crystal edge, with the "length" <0.1 that of crystal. The respec- 
tive scattering and losses over O.IL are much smaller. In this case the protons 
retain better chances for extraction with later passes than in the former case 
(perfect alignment). The secondary b of the scattered protons are still sufficiently 
large (~30 fim ^x'L here), so the "gap" x'L is not dangerous. 

We come to conclusion that a peak efficiency with imperfect crystal is achieved 
at some tilt x' ^0. In the real experiment one scans x' while searching the peak, 
and comes to this case automatically. We used the case brfiax/x^Li — 0.1 as an 
illustration; the optimal x' will be found automatically in the scan. Further 
on, we refer to this case as to the "pre-scatter" case, when protons first gently 
pre-scatter in the crystal edge to come later with low divergence but high b. 
Understandably, with imperfect crystal the prescatter case may appear also for a 
small negative tilt, x' <0. Then, F{x') may have two peaks, with a dip at x'=Q. 
The width of the dip at x' ~0 is also an indicator for bmax- 

3 Simulation procedure 

The crystal was located 61 m upstream of CO point of Tevatron lattice, with the 
edge at the horizontal distance of X=1.75 mm from the beam axis. At the crystal 
location, the machine parameters were /3a;=105.7 m, a^=0.109 (horizontally), 
and I3y=2l.b m, ay=Q.lAS (vertically); tunes (5a:=20.5853 and Qy=2Q.h7AA. The 
beam invariant rms emittance was 2.5 mm-mrad, which corresponds to vertical 
rms divergence 11.5 /irad and width 0.24 mm at the crystal location. 

Crystal was a Si(llO) slab 40 by 3 by 3 mm^, 0.64 mrad bent, with a perfect 
lattice, and curved with a constant curvature to deflect protons in vertical direc- 
tion. As an option, we model an amorphous layer at its edge and/or irregularities 
of the surface. The horizontal parameters ,t, x' of incident particles are defined by 
the mechanism of diffusion. The two processes, diffusion and crystal extraction, 
are unfold in E853. Beam parameters in the channeling plane (vertical) are not 
disturbed by the diffusion. The exact value of b^ax matters only w.r.t. t. Since 
t is unknown for the real crystal, we can postulate bmax=^ A^iHj and then model 
crystals with different t. 
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4 Results 



Fig. 2 shows the F{y') angular scan for x'=0. The peak F of an ideal crystal is 
~44 %. The same fig. shows scans for the crystals with t=l //m (i.e. t=hmax) a-^d 
i=50 /xm, where at y'=x'=0 F« 36 % and 32 % respectively. However, for the 
imperfect crystals the real peak was found at x' ^0 (Fig. 3). With optimized x', 
one has peak F ~ 42 % and 35 % for t=l fim and t=50 /um respectively. Notably, 
the efficiencies and scans arc quite weakly dependent on the crystal perfection. 
The width of y' scan was 50-55 /xrad fwhm in these cases. 




, , 20 
y (/xrad) 

Figure 2: F(y') scan for x'=0. Ideal crystal: (o) first-pass and (•) overall efficiencies. 
Imperfect crystal: overall efficiency with t—1 /um (★) and t—50 /im (*). 



The angular scan F{x') is in Fig. 3. The depth of the dip at x' ^0 (i.e. at 
perfect alignment) is ~14 % and ~7 % w.r.t. the peak for t=l jum and t=50 fim 
respectively. The width Ax' of the peculiarity (either peak or dip) near x' ~0 is 
roughly hmax/F which is 25 /Ltrad in our simulation. F is 1/2 of the maximum at 




Figure 3: F{x') scan near the peak; see fig.2 caption. 



x' ~14 mrad and —0.3 mrad for an ideal crystal {fwhm of the horizontal scan is 

~14 mrad), at x' ~15 mrad and -1.2 mrad for the crystal with t=l fim (/u;/im~16 
mrad), and at x' ~18 mrad and -5 mrad for the crystal with i=50 fj,m. (/u;/im~23 
mrad) . 
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The asymmetry of the scan, F{x') 7^ F{—x'), is due to the loss of the protons 
trapped in channeling near the crystal edge. With an ideal crystal, the asymmetry 
exists for any x' . With a septum width t, the asymmetry can be seen for an angling 
zizx' larger than t/L only. In our simulation with t=50 /7,m, the scan is symmetric 
indeed within ±1.3 mrad but asymmetric outside this range of x'; note that 50 
/xm/40 mm =1.25 mrad. We expect therefore this x'-threshold for an asymmetry 
to be a good measure of the septum width t. If one plots the magnitude of 
asymmetry, F{x')-F{—x'), as a function of x'L, he obtains a rough estimate of 
the beam distribution over b at crystal. The minimal step 6b = 6x'L=0.1 /im is 
much finer than the precision of coordinate detectors RiO.l mm ! 

Notice an abrupt decrease in F of the ideal crystal over the range of x'L from 
to —bmax- from 44 % at x'=0 to 28 % at x'= —bmax/L. This drop is an excellent 
opportunity to measure the primary bmax with a precision of (56=0.1 jim.. With 
an ideal crystal one can measure a distribution over the primary b (in the range of 
~1 /xm). With imperfect crystal in the same way one measures the distribution 
over secondary b (in a broad range from ~t to ~1 mm). Finally, the dependence 
F{x') for x'>0 gives actually the dependence of F on the septum width t ~ x'L. 

The distribution of the extracted particles over x at the crystal face is essential 
for understanding both the crystal interplay with other accelerator elements and 
the requirements for the crystal face perfection. We have found that one half 
of the extracted protons have penetrated into the crystal depth by >0.3 mm; 
another half had 6<0.3 mm. The y' divergence of the extracted beam was defined 
by the channeling properties of Si(llO) crystal; its full width 26c was Ril2.8 ^urad 
[0c— 6. 4: firad), and fwhm^9 /7,rad. The x' divergence was ~5 ^urad fwhm with the 
ideal crystal and ~12 /xrad fwhm with t=l /xm. It was increased due to scattering 
in inefficient passes. After the doublet of quadrupoles and the Lambertson-type 
magnet, two detectors (hodoscopes with 0.1 mm bins) were placed at 80.5 m (Dl) 
and 120.5 m (D2) downstream of the crystal to measure the bent-beam profiles. 
The horizontal profiles had width ~ 0.3 and 0.4 mm fwhm for the ideal and t=l 
yum crystals at Dl, and ~ 0.5 and 0.7-0.9 mm fwhm at D2. 

5 Optimization 

The extraction efficiency F is defined by the processes of channeling, scattering, 
and nuclear interaction in crystal. All the processes depend essentially on the 
crystal length L. Fig. 4 shows F{L). F is maximal, near 70 %, in L range 
0.4-1.0 cm, irrespective of the crystal perfection. 

6 Conclusions 

The key information of the multi-pass crystal extraction can be obtained from 
the analysis of the horizontal angular scan of efficiency. In the considered way 
one can study the impact parameters of halo particles and the structure of the 
crystal edge with an accuracy as fine as 0.1 ^m. 

The extraction efficiency is expected as high as ~40 % irrespective of the 
crystal septum width, and can be increased up to ~70 % with the use of a shorter 
(<1 cm) crystal. The difference in efficiency between the ideal and imperfect 



5 



F- 

8 o 

0.6 - 

o 

• o 

0.4 - 

9 • 

0.2 - 

: 1 2 3 4 
L (cm) 

Figure 4: F{L) for ideal (o) and t=l //m (•) crystals. 



crystals is very low, because of predominance of the multi-passes in extraction at 
high energies, and partly due to the found effect of a gentle " prescattering" in 
the edge of a crystal tilted horizontally. This provides an elementary solution to 
the problem of a finite septum width and infinitesimal impact parameters. 

One general trend in the results of simulations, from SPS [7] to Tevatron to 
LHC [6] , is worthwhile to mention: the difference in efficiency of the ideal crystal 
and crystal with imperfect surface vanishes with energy E, because the scattering 
angle reduces faster than 9c does 
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